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ABSTRACT 


A numerical method for computing three-dimensional, time-dependent 
incompressible flows is presented. The method is based on a fractional- 
step, or time-splitting, scheme in conjunction with the approximate- 
factorization technique. It is shown that the use of velocity boundary 
conditions for the intermediate velocity field can lead to inconsistent 
numerical solutions. Appropriate boundary conditions for the interme- 
diate velocity field are derived and tested. Numerical solutions for 
flows inside a driven cavity and over a backward-facing step are pre- 
sented and compared with experimental data and other numerical results. 
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1 . INTRODUCTION 


In this paper we present a numerical method for solving three- 
dimensional, time-dependent incompressible Navier-Stokes equations. The 
major difficulty in obtaining a time-accurate solution for an incompres- 
sible flow arises from the fact that the continuity equation does not 
contain the time-derivative explicitly. The constraint of mass conser- 
vation is achieved by an implicit coupling between the continuity equa- 
tion and the pressure in the momentum equations. One can use an explicit 
time-advancement scheme which obtains the pressure at the current time- 
step such that the continuity equation at the next step is satisfied. 
However, for fully implicit or semi-implicit schemes, the aforementioned 
difficulty prevents the use of the conventional alternating-direction- 
implicit (ADI) scheme to advance in time as is the case for compressible 
flows. This difficulty can be avoided in two-dimensional cases by 
reformulating the problem in terms of the vorticity and stream-function. 
For three-dimensional problems, one can introduce an artificial compres- 
sibility into the continuity equation to include the required time- 
derivative for an ADI scheme. This is satisfactory, however, only for 
the steady— state solutions [1]. For unsteady problems, since the effect 
of the artificial compressibility has to be minimized, this approach 
produces a highly stiff system for numerical solutions [2]. 

The objective of the present work is to develop a numerical method 
for solving the incompressible Navier-Stokes equations satisfying the 
conservation of mass exactly (within machine round-off). It will be 
also required that the numerical scheme preserve the global conservation 
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of momentum, kinetic energy, and circulation in the absence of time- 
differencing errors and viscosity. It can be shown that failure to 
preserve these conservation properties can lead to numerical instabili- 
ties [3], In order to stabilize the calculations using methods that do 
not preserve these properties, artificial viscosity is often introduced 
either explicitly or implicitly by using dissipative finite-difference 
schemes, especially for high-Reynolds-number flows. 

The method developed herein is based on a fractional-step method 
(e.g., [4,5]) in conjunction with the approximate-factorization tech- 

nique [6,7]. The flow field is represented on a staggered grid [8]. 

The problem of concocting boundary conditions for the intermediate 
(split) velocity field is addressed and it is shown that the use of 
velocity boundary conditions can lead to inconsistent and erroneous 
results. Appropriate boundary conditions for the intermediate— velocity 
field are derived using a technique similar to that of LeVeque and 
Oliger [9]. The Poisson equation for the pressure correction is solved 
by a direct method based on trigonometric expansions. In this way the 
continuity equation is satisfied to machine accuracy at every time-step. 

The numerical procedures used in the present method are described 
in Section 2. Section 3 provides a derivation of the boundary condi- 
tions for the intermediate-velocity field, and numerical results for two 
different flow geometries are presented in Section 4; a summary is given 
in Section 5. 
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2. NUMERICAL METHOD 


In this section we present a variant of the fractional-step method 
used by Chorin [4] for time-advancement of the Navier-Stokes and continu- 
ity equations for incompressible viscous flows: 
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3x. 

3 


Vj 


Js- + _L JL-l- u. , 
3x ± Re 3x^ 3Xj l 



(1) 

( 2 ) 


Here, all variables are nondimensionalized by a characteristic velocity 
and length scale, and Re is the Reynolds number. 

The fractional step, or time-split method, is in general a method 
of approximation of the evolution equations based on decomposition of 
the operators they contain. In application of this method to the Navier- 
Stokes equations , one can interpret the role of pressure in the momentum 


equations as a projection operator which projects an arbitrary vector 
field into a divergence-free vector field. A two-step time-advancement 
scheme for Eqs . (1) and (2) can be written as follows: 

G i " U i 1 nH n . H n-i } + I _L (1L + *L + *L\ (u. 4- u n ) , (3) 

~ 2 ( i H ‘ ’ 6 1 1 


n+i 


U. 


- U 


At 


— = -G((f> n+1 ) 


W) 


with 

D(uf ■) - 0 , <5) 

where H i = -(S/fix^u^ is the convective terms, f is a scalar to be 
determined, <5/6x i represents discrete finite difference operators, 
and G and D represent discrete gradient and divergence 
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operators, respectively. We used the second-order-explicit Adams- 
Bashforth scheme for the convective terms and the second-order-implicit 
Crank-Nicolson for the viscous terms. Implicit treatment of the viscous 
terms eliminates the numerical viscous stability restriction. Equa- 
tion (3) is a second-order-accurate approximation of Eq. (1) with Sp/Sx^ 
excluded. By substituting Eq. (4) into (3), one can easily show that 
the overall accuracy of the splitting method is still second order. Note 
that 4> is different from the original pressure: in fact, 
p = <f> + (A t/2Re) V 2 <j> . All the spatial derivatives are approximated with 
second-order central differences on a staggered grid [8], Figure 1 
illustrates the staggered grid. With the staggered mesh, the momentum 
equations are evaluated at the velocity nodes, and the continuity equa- 
tion is enforced for each cell. One important advantage of using the 
staggered mesh for incompressible flows is that ad hoc pressure boundary 
conditions are not required. Furthermore, it can be shown [10] that 
with this approximation of spatial derivatives and in the absence of 
time-differencing errors and viscosity, global conservations of momen- 
tum, kinetic energy, and circulation are preserved. 

Equation (3) can be rewritten as 

(1 - A 1 - A 2 - A 3 ) (u^ - u*) = -y (3 hV - h"' 1 ) 

+ 2 (A x + A 2 + A 3 )u” , (6) 

where A 1 = (At/2Re) (6 2 /5x 2 ) , A 2 = (At/2Re) (6 2 /6x 2 ) , A 3 = (At/2Re) (6 2 /<5x 2 ) . 
The left-hand side of Eq. (6) is then approximated as follows: 
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( 7 ) 


(1 - Aj) (1 - AjHl - A 3 ) (u ± - u£) - ^(3 hJ - h" _1 ) 

+ 2(A X + A 2 + A 3 )u" . 

Equation (7) is an 0(At 3 ) approximation to Eq. (6). However, it requires 
inversions of tridiagonal matrices rather than inversion of a large 
sparse matrix, as in the case of Eq. (6). This results in a significant 
reduction in computing cost and memory. 

Equations (4) and (5) can be solved as coupled system equations 
for u n+1 and (}> n+1 with boundary conditions for u? +1 . Note that since 


<f> 


n+i 


is defined at the center of each cell, there is a sufficient num- 


ber of equations for u* +1 and <J> n+1 without the need for boundary con- 


dition for (J) 
n+i 


n+i 


Equations (4) and (5) can be combined to eliminate 

, n*f i 


u and thus obtain a set of equations for . For the cells not 

i 

adjacent to the boundaries, these equations take the form of the discrete 
Poisson equation. 


(— 

\6xi 


6 2 

2 2 
6x 2 6x 3 


) 


t> n+1 (i,j >k) = Du 


= Q(i , j ,k) , (8) 

for i = 2,3,...,!^ - 1, j = 2,3 N 2 - 1, k = 2,3,. ...N 3 - 1. For 

the cells adjacent to the boundaries, incorporation of the velocity 
boundary conditions yields a modified set of equations. For example, for 
the cells adjacent to the lower boundary (j =1), we obtain 


\5x 2 <5x 3 / 


$ n+1 (i,l,k) + 


1 


x 2 (f) " x 2 (y) 


<t n+1 (i,2,k) - <t. n+1 (i,l,k) 
x, (2) - x 2 (l) 


1 A 1 u 2 +1 ( 1, 2 ,k ) U2 ( x * 2 ,k ) 

At Du " At 


X 2(l) " X 4l) 


s Q(i.l.k) 


(9) 
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for i = 2,3 N x - 1, k = 2,3,...,N 3 - 1. A solution to Eq. (8) and 

the corresponding boundary equations can be easily obtained using trans- 
form methods [11]. Let 

Nj-i N 3 -i 

$ n+1 (i>j»k) = $(«.,j ,m)cos 

£=o m=o 

( 10 ) 



for i = l,2,...,N lf j = 1,2,...,N 2 , k = Here we have 

assumed that uniform mesh spacings are used in the streamwise and span- 
wise directions (x 1 and x 3 ). If nonuniformly spaced mesh points are 
used in these directions, multigrid methods appear to be the best 
alternative to the transform method used here. Substituting Eq. (10) 
and the corresponding expansion for Q into (8) and (9) and using the 
orthogonality property of cosines, we obtain 


^ - k^ - k%j> = Q(£, j ,m) , (11) 

6x 2 

where k^ = 2 [ 1 - cos (ttJI/N x ) ] /Ax* and k^ = 2 [ 1 - cos (7rm/N 3 ) ] / Ax* are 
the modified wave numbers. For each set of wave numbers, the above tri- 
diagonal system of equations can be easily inverted, and 4> n+ is 
obtained from Eq. (10). The final velocity field u" +1 is then obtained 
from 


n+i 

u. 


u ± - AtG(<t» n+1 ) 


( 12 ) 


3. BOUNDARY CONDITIONS 

Boundary conditions for the intermediate velocity fields in time- 
splitting methods are generally a source of ambiguity. At each complete 
time-step, only the boundary conditions for the velocity field are given 
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and those of the intermediate velocity field are unknown. We will show 
in this section that except when the boundary conditions for the inter- 
mediate velocity field are chosen to be consistent with the governing 
equations, the solution may suffer from appreciable numerical errors. 

In the present work, we derive the appropriate boundary conditions for 
the intermediate velocity field using a method suggested by LeVeque 
and Oliger [9]. 

To construct the proper boundary conditions for u^, we regard (k 
as an approximation to u?(x,t n+i ) where the continuous function 
u*(x,t) satisfies 

3u 


i _ , 1 8 9_ * 

at H i Re ax. 3x. u i 
J J 




( 13 ) 


Hence 


u. = u.(x,t + At) 
l l - n 


3u: 


a 2 u* 

= u«(x,t n ) + At-g^ + y At 2 —± + .., 

9 1 


- + 4t ( H i + Te V2u ?) + I i ( H J + Te ^tj 


Since u?(x,t ) = u. (x,t ), 
i - n i - n 


“i ■ u l ( S- t n ) + 4t ( H i + R? ^“i) + T “ 2 Ji ( H i + ite v2u i) + " 

(-5T + 1^) + i “ 2 (-r + £ 4 p) - - 


= u. (x,t ) + At 
i - n 


(14) 


( 15 ) 


= u i ( ^’ t n+i ) + At ^- + Q(At 2 ) . 
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By keeping the first two terms, we have boundary conditions accurate to 
Q(At 2 ). Since p = 4> + 0(At/Re) , we can in fact use = u? +1 + At 


dx. 

X 


with the same accuracy, thus avoiding the computation of pressure 
explicitly. 

These boundary conditions are tested in computing the following 
two-dimensional unsteady flow which is a solution to the Navier-Stokes 
and continuity equations [12,13]: 

u 1 (x 1 ,x 2 ,t) = -cos x 1 sin x 2 e 

— 2 t 


— 2 t 


u 2 (x 1 ,x 2 ,t) = sin x 1 cos x 2 e 

1 - 4 1 

p(x 1 ,x 2 ,t) = - (cos 2x 1 + cos 2x 2 )e 


(16) 


The maximum difference between the exact and numerical solutions from 
four different runs are listed in Table I. The superiority of the results 
using boundary conditions (15) over the results using the velocity bound- 
ary condition, u^ = u? +1 , is clearly evident. In fact, the error for 
the latter case increases when the mesh is refined. This indicates that 
the fractional-step method with the latter boundary condition is an 
inconsistent scheme. To determine the overall accuracy of the scheme 
using the boundary condition (15), three calculations are performed 
with three different mesh sizes but keeping the Courant number constant. 
The variation of the maximum error with mesh refinement is plotted in 
Fig. 2; it shows that the scheme is indeed second-order accurate. 


4 . NUMERICAL EXAMPLES 


In this section we present numerical results obtained from applica- 
tions of the aforementioned numerical method to two laminar flow problems. 
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Both problems have been used widely as standard test cases for evaluating 
the stability and accuracy of numerical methods for incompressible flow 
problems. 

4.1 Flow in a Driven Cavity 

Figure 3 shows the geometry and the boundary conditions for the 
flow in a driven cavity together with the appropriate nomenclature. 

Flow is driven by the upper wall, and several standing vortices exist 
inside the cavity whose characteristics are functions of Reynolds numbers. 
Figures 4 and 5 show the computed results of streamlines, contours of 
constant vorticity, and velocity vectors for several Reynolds numbers. 

The purpose of the velocity-vector figures is to show the corner eddies , 
which are too weak to be displayed clearly by the streamlines. They are 
drawn parallel to the flow direction at each mesh point. At Re = 1 , 
this flow is almost symmetric with respect to the centerline, and two 
corner eddies are visible. As Reynolds number increases, the center of 
the main vortex moves toward the downstream corner before it returns 
toward the center at higher Reynolds number. In the Reynolds number 
range, 1000-2000, the third corner eddy is formed at the upper left 
corner. At Re = 5000, a tertiary corner eddy is visible. In Fig. 6, 
the velocity at the middle of the cavity for Re = 400 is shown in 
comparison with other computed results. Two numerical results with 
different grid sizes are shown from the present computations. Both 
results, 21 x 21 and 31 x 31, are in good agreement with those of 
Burggraf [14] and Goda [15]. Although not shown here, Goda reported a 
rather poor agreement when he used a 21 x 21 mesh for this Reynolds 
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number. In Table II, the magnitudes of the stream-function and vorticity 
at the center of the primary vortex from the present calculations are 
compared with those of the other investigators [16-18] at Re = 1000. 

In Table III, for different Reynolds numbers, the same quantities from 
the vorticity stream-function calculations of Ghia et al. [16] and 
Schreiber and Keller [18] are compared. 

In their experimental study, Koseff et al. [19] observed Taylor- 
Gortler type vortices, which are formed as a result of the streamline 
curvature owing to the primary vortex. This is the first observation of 
such vortices in a cavity flow. Their numerical simulation, however, 
failed to reproduce this three-dimensional structure. To initialize the 
calculation in the present study, small random disturbances in the span- 
wise direction (x 3 ) were added to the solution of two-dimensional cases. 
Using periodic boundary conditions in the spanwise direction, the com- 
putations were carried out for various Reynolds numbers. In Fig. 7, the 
velocity vectors in the plane perpendicular to the primary vortex show 
the existence of the counterrotating vortices at Re = 1000. Although 
Goda [15] calculated the flow in a three-dimensional cavity, no such 
three-dimensional structure was reported. 

4.2 Flow over a Backward-Facing Step 

The flow over a backward-facing step in a channel provides an 
excellent test case for the accuracy of numerical method because of the 
dependence of the reattachment length x^ on the Reynolds number. 
Excessive numerical smoothing in favor of stability will result in fail- 
ure to predict the correct reattachment length. 
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The geometry and boundary conditions for this flow are shown in 
Fig. 8. At the inflow boundary, located at the step, a parabolic profile 
was prescribed. Both Neumann and Dirichlet outflow boundary conditions 
were used, and the two results are identical. In Fig. 9, numerical 
results for different Reynolds numbers are shown in comparison with the 
experimental and computational results of Armaly et al. [20]. The 
dependence of the reattachment length on Reynolds number is in good 
agreement with the experimental data up to about Re = 500. At Re = 600, 
the computed results start to deviate from the experimental values. A 
mesh-refinement study, as well as variation of the location of downstream 
boundary at this Reynolds number, showed that the difference between the 
experimental and computational results is not a result of numerical 
errors. It is most likely, as Armaly et al. [20] have pointed out, that 
the difference is due to the three-dimensionality of the experimental 
flow at this Reynolds number. In comparison with the numerical results 
of Armaly et al. [20] (using TEACH code), however, the present results 
show a much higher reattachment length. 

Armaly et al. [20] reported the existence of a secondary separation 
bubble on the no-step wall at Re = 1000. The length of the secondary 
bubble at Re = 1000 was 10.4 step-heights and the length decreased for 
higher Reynolds numbers. Figure 10 shows the computed streamlines at 
Re = 600, indicating the secondary separation bubble on the no-step wall; 
the bubble length is 7.8 step-heights. At Re = 800, the length increased 
to 11.5 step-heights. 
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SUMMARY 


A numerical method was presented for solving three-dimensional , 
time-dependent incompressible flows; the method is based on the fractional- 
step method used in conjunction with the approximate-factorization scheme. 
The three-dimensional Poisson equations was solved directly by a trans- 
form method, and the velocity field satisfied the continuity equation up 
to machine accuracy. The method is second-order accurate in both space 
and time. Proper boundary conditions for the intermediate (split) veloc- 
ity field were derived and tested against a known solution, and laminar 
flows in a driven cavity and over a backward-facing step were calculated 
at several Reynolds numbers. The numerical results are in good agree- 
ment with experimental data and other numerical solutions. 
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TABLE I 


Maximum Error after 30 Steps: e /u 

r max max 



Boundary conditions 

Grid 

- n+i 

u i = u i 

* n+i , A ^ 

u i ■ u i + ^ 9Xi 

o 

CM 

X 

O 

CM 

8.172 x 10 -4 

.1.085 x 10 -4 

40 x 40 

1.127 x io“ 3 

7.678 x 10 -5 


TABLE II 


Stream-Function and Vorticity at Center of Primary Vortex at Re = 1,000 



Ghia et al. [16] 

Present 

Benjamin and Denny [17] 

Schreiber and Keller [18] 


(129 x 129) 

(97 x 97) 

(101 X 101) 

(141 x 141) 

*c 

-0.118 

-0.116 

-0. 118 

-0.116 

? c 

-2.050 

-2.026 

-2.044 

-2.026 
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TABLE III 

Stream-Function and Vorticity at Center of Primary Vortices 
for Different Reynolds Numbers 


Re 

Present 

<«c> 

Ghia et al. [16] 
*c «c> 

Schreiber and Keller [18] 
*c «c> 

i 

-0.099 

(-3.316) 

■ 

— 

-0.100 

(-3.232) 


65 

x 65 



121 

x 121 

100 

-0.103 

(-3.177) 

-0.103 

(-3.166) 

-0. 103 

(-3. 182) 


65 

x 65 . 

129 

x 129 

121 

x 121 

400 

-0.112 

(-2.260) 

-0.114 

(-2.295) 

-0.113 

(-2.281) 


65 

x 65 

257 

x 257 

141 

x 141 

1,000 

-0. 116 

(-2.026) 

-0.118 

(-2.050) 

-0.116 

(-2.026) 


97 

x 97 

129 

x 129 

141 

x 141 

3,200 

-0. 115 

(-1.901) 

-0.120 

(-1.989) 

- 

— 


97 

x 97 

129 

x 129 



4,000 

-0.114 

(-1.879) 

- 

— 

-0.112 

(-1.805) 


97 

x 97 



161 

x 161 

5,000 

-0.112 

(-1.812) 

-0.119 

(-,1.860) 

- 



97 

x 97 

257 

x 257 
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Figure Captions 


FIG. 1. The staggered mesh in two dimensions. 

FIG. 2. Maximum error as a function of mesh refinement. 

FIG. 3. Geometry of the driven cavity flow. 

FIG. 4. Streamlines and contours of constant vorticity. (a) Re = 1; 
(b) Re = 400; (c) Re = 2,000. 

FIG. 5. Velocity vectors, (a) Re = 1; (b) Re = 100; (c) Re = 400; 
(d) Re = 1,000; (e) Re = 2,000; (f) Re = 5,000. 

FIG. 6. Profile of streamwise velocity at the midplane of the cavity 
(x x = 0.5L) for Re = 400. 

FIG. 7. Velocity vectors in an (x 2 - x 3 ) plane through the geometric 
center of the cubic cavity. 

FIG. 8. Flow over a backward-facing step (1:2 expansion ratio). 

FIG. 9. Reattachment length as a function of Reynolds number. 

FIG. 10. Streamlines at Re = 600. 


20 




21 












Fig. 


(e)(f) 





>7V“t- 

-t- 

T 

T7 

— 

— 


X 

T-t- 

l 

7 7 ; ^ — 


1 1 1 

1 

1 

/ 

/ 

/ 

/ /\ X \ 

\ 

\ 1 

1 1 

I- 

1 / 

/ 

/ 

/ 

/ /xX \ \ \ 

( 1 1 

1 

1 

\ 

\ 

\ 

\ 

1 

1 

l 1 

f 

1 t 

1 1 

\ 

\ \ 

\ 

\ 

\ 

\ 

1 

1 

1 1 i 1 

/ 1 1 

1 


\ 

\ 

\ 

\ 


1 

/ / 

i 

/ / 

1 1 

\ 

\ \ 

\ 

\ 

\ 

\ 

\ 

1 

1 / / / 

/ / / 

/ 

[ 

\ 

1 

\ 

1 

\ 


1 

/ / 

/ 

1 1 

1 \ 

\ 

\ \ 

\ 

\ 

\ 

\ 

\ 

1 

I / / / 

/ / / 
1 1 1 

/ 

/ 

1 

1 


1 

1 1 

1 

1 1 

\ \ 

\ 

\ \ 

\ 

\ 

\ 

\ 

\ 

1 

till 

/ 

/ 

f 

/ 

/ 

1 


1 

1 1 

\ 

\ \ 

\ 1 

1 

1 1 

1 

1 

l 


i 

1 

Mil 

1 1 1 

/ 

1 

/ 

/ 

1 

1 


1 

\ \ 

\ 

\ \ 

\ 1 

1 

1 1 

1 

1 

1 

1 

1 

1 

M. \ \ 
M \ \ 


1 

1 

/ 

/ 

1 

1 


1 

\ \ 

\ 

\ \ 

\ 1 

1 

1 / 

1 

t 

1 

f 

f 

1 

1 

1 1 1 

1 

1 

/ 

/ 

/ 

1 

1 

1 


1 

\ \ 

\ 

\ \ 

\ 1 

1 

1 I 

1 

1 

/ 

1 

1 

\ \ \ \ 

\ 1 | 

| 

\ 

/ 

( 


1 

1 \ 
1 1 

\ 

\ 1 

1 1 

1 

1 

1 

t 

1 

/ 

1 

1 

M \ \ 

\ | [ 

1 

| 

/ 

| 

1 

l 


1 

| 

1 1 

! 1 

1 

i 1 

1 

1 

1 

1 

1 

1 

\ \ \ \ 


1 

I 

1 

| 

1 

| 

1 

| 

1 

| 

1 


1 

1 

1 1 
1 1 

1 

I 

1 1 

1 1 

l 1 
1 1 

1 

1 

1 1 
1 1 

1 

1 

l 

1 

1 

1 

I 

1 

1 

1 

1 

1 

\ \ \ \ 
MM 

1 I 1 
1 1 1 

1 

| 

1 

1 

1 

| 

1 

1 

1 

\ 

1 

1 


/ 

1 

1 1 
1 i 

1 

1 

1 1 
1 1 

1 1 
1 1 

1 

1 

1 1 
1 1 

1 

\ 

1 

\ 

1 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

Mil 

MM 

/ 1 1 

1 

\ 

\ 

\ 

\ 

1 


1 

/ / 

/ 

1 l 

f 1 

1 

1 \ 

\ 

\ 

\ \ \ 

\ 

1 III 

/ / ! 

\ 

\ 

\ 

\ 

\ 

\ 


1 

/ / / / / 

/ 1 

\ 

\ \ W\\ \ 

\ 

f / / / 

// / 

\ 

\ \ \ 

\ 

\ 


/ 

/////! \ 

\ 

\ X 




X 

\ 

\ 

i / / ^ 


x 2 

4 


w\ / / 

\ \ ' i 
\ \ i / 

\ \ i / //" 

mi///" 

\ i i i / 
iiiiii/ 
/ ii \ \ \ 

K / \ \ 


\ \ 
\ \ 
\ \ 
\ \ 
\ \ 
\ \ 
i i 

L // 


/ / / /_. 
/ ( / 
i / 
i / 
i / 
i / 


i / 





\ \ i 
N \ \ i 
"W \ i 
\ \ i 
\\ \ i 1 
i MM 

v J ^ / / / / 

s / / 
^ ^ / 


/ s'- 
/ //^ 

/ / 

/ / // 

I I / / 

/ /./ / 

I I I I / 

I \ \ \ ■ 
\ \ v 
X ^ 


'S 


.x \ \ / / 

_\ \ W / 

_ \ \ l / / 

_x \ W / 

.x \ l / / 

^ \ \ I / / J. 
Mil//. 

_ / / ( \ \ X ^ 

-/ / \x-M 

-MXx 


x 3 


Re = 1000 


Fig. 7 


29 



U) 

o 



1C INFLOW u-, = U 2 = 0 ON WALLS 

'^Z/77Z7/A 


DIVIDING STREAMLINE 


DOWNSTREAM 

BOUNDARY 

! 3 


REATTACHMENT 


a^ u i 


i ox i 


! 3 


777777777 :. 


F 


x,/h 



Fig. 9 


31 





1. Report No. 

NASA TM-85898 


2. Government Accession No. 


3. Recipient's Catalog No. 


4. Title and Subtitle 

APPLICATION OF A FRACTIONAL-STEP METHOD TO 
INCOMPRESSIBLE NAVIER-STOKES EQUATIONS 


5. Report Date 

March 1984 


6. Performing Organization Code 

ATP 


7. Author(s) 

J. Kim and P. Mo in 


9. Performing Organization Name and Address 

Ames Research Center 
Moffett Field, CA 94035 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, DC 20546 


8. Performing Organization Report No. 

A-9665 


10. Work Unit No. 

T-6465 


11. Contract or Grant No. 


13. Type of Report and Period Covered 

Technical Memorandum 


14. Sponsoring Agency Code 

505-31-01-01 


15. Supplementary Notes 


Point of Contact: J. Kim, Ames Research Center, MS 292A-1, Moffett Field, 

CA 94035 (415) 965-5576 or FTS 448-5576 

16. Abstract 

A numerical method for computing three-dimensional, time-dependent 
incompressible flows is presented. The method is based on a fractional- 
step, or time-splitting, scheme in conjunction with the approximate- 
factorization technique. It is shown that the use of velocity boundary 
conditions for the intermediate velocity field can lead to inconsistent 
numerical solutions. Appropriate boundary conditions for the intermediate 
velocity field are derived and tested. Numerical solutions for flows 
inside a driven cavity and over a backward-facing step are presented and 
compared with experimental data and other numerical results. 


17. Key Words (Suggested by Author(s)) 

Numerical method 
Navier-Stokes equations 
Fractional-step methods 


18. Distribution Statement 

Unlimited 

Subject category: 64 

19. Security Classif. (of this report) 

20. Security Classif. (of this page) 

21. No. of Pages 

22. Price* 

Unclassified 

Unclassified 

35 

A0 3 


*For sale by the National Technical Information Service, Springfield, Virginia 221 61 






















